PENTAGONAL DOMAIN EXCHANGE 



SHIGEKI AKIYAMA AND EDMUND HARRISS 



Abstract. Self-inducing structure of pentagonal piecewise isometry is 
applied to show detailed description of periodic and aperiodic orbits, 
and further dynamical properties. A Pisot number appears as a scaling 
constant and plays a crucial role in the proof. Further generalization is 
discussed in the last section. 

Adler-Kitchens-Tresser [1] and Goetz [10] initiated the study of piecewise 
isometries. This class of maps shows the way to possible generalizations of 
results on interval exchanges to higher dimensions 130] , In this paper 
we examine the detailed properties of the map shown in Figure [1] from an 
algebraic point of view. 




Figure 1. A piecewise rotation T on two pieces. The tri- 
angle is rotated 2tt/5 around a and the trapezium is ro- 
tated 27r/5 around b. Periodic points with short periods are 
shown below, in two colours to illustrate that they cluster 
into groups, each forming a pentagon. 



The goal of this paper is to see how this map is applied to show number 
theoretical results. First we reprove that almost all orbits in the sense of 
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Lebesgue measure are periodic, and in addition, there are explicit aperiodic 
points. Second we show that aperiodic points forms a proper dense subset of 
an attractor of some iterated function system and are recognized by a Biichi 
automaton (c.f. Figure [14]) . The dynamics acting on this set of aperiodic 
points are conjugate to the 2-adic odometer (addition of one) whose exphcit 
construction is given (Theorem [3]) . As a result, we easily see that all aperi- 
odic orbits are dense and uniformly distributed in the attractor. We finally 
give a characterization of points which have purely periodic multiplicative 
coding by constructing its natural extension (Theorem [6]). In doing so we 
obtain an intriguing picture Figure [16] that emerges naturally from taking 
algebraic conjugates, whose structure is worthy of further study. We discuss 
possible generalizations for 7-fold and 9-fold piecewise rotations in Section 

m 

A dynamical system is self-inducing if the first return map to some sub- 
set has the same dynamics as the full map. The most important example is 
the irrational rotation, presented as exchange of two intervals. An elemen- 
tary example begins with $ shown in Figure [2] For this interval exchange. 



I 



(D(B) = B - (D(A) = A * I - A-i 

Figure 2. An interval exchange map ^, where A = 



now consider the second interval B. As shown in Figure [3] this interval is 
translated to the left once, and to the right. Thus ^^(Bi) is back in B, the 
interval B2 requires one more step, but <I>^(i?2) also lies within B. This first 
return dynamics on B is therefore conjugate to the dynamics on AUB. Self- 
inducing subsystem of two interval exchange corresponds to purely periodic 
orbits of continued fraction expansion and they are efficiently captured by 
the continued fraction algorithm. 

This gives a motivation to study the interval exchange transform (lET) 
of three or more pieces, trying to find higher dimensional continued fraction 
with good Diophantine approximation properties. The study of self- inducing 
structure of lET's was started by a pioneer work of Rauzy [26|, now called 
Rauzy induction, and got extended in a great deal by many authors including 
Veech [3l] and Zorich p3] . see |32j for historical developments. 

Self-inducing piecewise isometries emerged from dynamical systems as a 
natural generahzation of lET [SH [T] [5] [IH [IS] [H] [22] and the first return 
dynamics appears in outer billiards [28l[6]. Like lET they provide a simple 
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Figure 3. The interval exchange $ is self-inducing. The 
intervals Bi and B2 are swapped by the first return map of 
$ on the interval B. 

setting to study many of the deep and perplexing behaviors that can emerge 
from a dynamical system. 

The self-inducing structure links such dynamical systems to number the- 
oretical algorithms, such as, digital expansions and Diophantine approxima- 
tion algorithms, and allows us to study their periodic orbits by constructing 
their natural extensions. This idea leads to complex and beautiful fractal 
behavior. 

Our target is the piecewise isometry in Figure [H but to illustrate the 
bridge formed between the two fields let us begin with a simple conjecture 
from number theory: 

Conjecture 1. For any —2 < A < 2, each integer sequence defined by 
< ttn+i + Aa„ -|- a„„i < 1 is periodic. 

Since a„+2 G Z is uniquely determined by (a„,a„4-i) G Z^, we treat this 
recurrence as a map (a„,an+i) 1— )■ (an+i,a„+2) acting on Z^. It is natural 
to set A = — 2cos(0) to view this map as a 'discretized rotation': 



with eigenvalues exp(ziz\/—19). As the matrix is conjugate to the planar 





1 
cos 6 — sin 
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where {x) is the fractional part of x. Therefore this gives a rotation map 
of angle 9 acting on a lattice PZ^ but the image requires a bounded per- 
turbation of modulus less than two to fit into lattice points of PZ^. For 
conjecture [U we expect that such perturbations do not cumulate and the 
orbits stay bounded, equivalently, all orbits become periodic. 

A nice feature of the map (a„,a„+i) i— )• (a„+i,a„4.2) is that it is clearly 
bijective on 7? by symmetry, while under the usual round off scheme, the 
digital information should be more or less lost by the irrational rotation. 
This motivates dynamical study of global stability of this algorithm. 

The conjecture is trivial when A = 0, ±1. Among non-trivial cases, 
the second tractable case is when 6 is rational and A is quadratic over Q. 
Akiyama, Brunette, Petho and Steiner [3] proved: 

Theorem 1. The conjecture is valid for A = ^^^^ , ±\/2, 

It seems hard to prove Conjecture [1] for other values. The case A = 
was firstly shown by Lowenstein, Hatjispyros and Vivaldi [2T] with 

heavy computer assistance. A number theoretical proof for appeared 
in [2j, whose proof is short but not so easy to generalize. We try to give an 
accessible account using self-inducing piecewise isometry in the case X = uj = 
, together with its further dynamical behavior. The proof in Section [T] 
is basically in [3]. However this version may elucidate the background idea 
and is directly connected to the scaling constant of self-inducing structure 
of piecewise isometry acting on a lozenge. 

A Pisot number is an algebraic integer > 1 whose conjugates have modulus 
less than 1. Throughout the paper, we will see the importance of the fact 
that the scaling constant of self-inducing system is a Pisot number. Our all 
discussions heavily depend on this fact. Indeed, Pisot scaling constants often 
appear in self-inducing structures of several important dynamical systems, 
for e.g., lET and substitutive dynamical systems. We discuss this point 
in Section O It is pretty surprising that we see this phenomenon in cubic 
piecewise rotations as well. We hope this paper gives an easy way to access 
this interesting area of mathematics. 

We wish to show our gratitude to P.Hubert, W.Steiner and F.Vivaldi 
for helpful comments and relevant literatures in the development of this 
manuscript. 



1. Proof of the periodicity for golden mean 

Setting C = exp(27ri/5), we have uj = — — ^-nd l/co = ( + The 
integer ring of Q(C) coincides with the ring generated by C in Z, Z[(^] is 
a free Z-module generated by 1,C,C^)C^- Hereafter we use a different base 
as a Z-module: 

Lemma 2. Z[("] is a free Z-module of rank 4 generated by l,a;, C^^C- 
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Proof. From uj = —(^'^ — C ^, we have 

xi + X2UJ + {yi + y2Uj)C = {xi + ^2) + [vi + 2/2)C + (2/2 - X2)C'^ - X2C^. 
On the other hand 

ao + oiC + a2C^ + osC^ = (oo - 02 + 03) - ^3^ + ((ai - 02 + 03) + (02 - a3)uj)(. 

□ 

Taking the complex conjugate, the same statement is vaHd with another 
basis l,uj, ,uj(~^ . Thus each element in has a unique expression: 

x-C""^y {x,y£Z[uj]). 

Denote by (x) the fractional part of x G M. Then a small computation gives 

< On + a;a„+i + a„+2 < 1 

On + oja-n+i + a„+2 = {uan+i) 

(uan) - -{u}an+i) + {ujan+2) = (mod Z) 

Xn- iC + C~^)Xn+l + Xn+2 = (mod Z) 

(x„+i - C~"^2;„+2) = C^ixn - C^Xn+i) (mod C^Z) 

and x„ = (a;a„). Our problem is therefore embedded into a piecewise isom- 
etry T acting on a lozenge [0, 1) + (— C~^)[0, 1): 




x/C Im(x/C) > 

(x-l)/C Im(x/C)<0' 



The action of T is geometrically described in Figure [TJ The lozenge L = 
[0, 1) + (— C~^)[0, 1) is rotated by the multiplication of — and then the 
trapezoid Z which falls outside L is pulled back in by adding — C~^- In 
total, the isosceles triangle A is rotated clockwise by the angle Svr/S around 
the origin and the trapezoid Z is rotated by the same angle but around 
the point ^ + i V5(5+2v^ ^ q ^ _|_ Q_ggg2f indicated by a black spot, that 
is the intersection of two diagonals. Our aim is to show that each point 
x G n L gives a periodic T-orbit. 

A. Goetz [TU] gave a slightly different map. Ours is an 'inclined' modifi- 
cation of p!6] and [3]. 

Clearly the map T is bijective and preserves 2-dimensional Lebesgue mea- 
sure fi. However the measure dynamical system {L,iy,M,T) (with the a- 
algebra B of Lebesgue measurable sets) is far from ergodic. It turned out 
that orbits of T is periodic for almost all points but for an exceptional set 
of Lebesgue measure zero. Our goal is to prove that the set Z[C] has no 
intersection with this exceptional set. This is not so obvious since Z[C] is 
dense in L because Z[a;] is dense in R. 

To illustrate the situation, it is instructive to describe an orbit of 1/3. 
See Figure HI 
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Figure 4. The orbit of 1/3 



Later we wih show that the orbit of 1/3 is aperiodic and forms a dense 
subset of the exceptional set of aperiodic points. Roughly speaking, our task 
is to show that Z[(^] n L has no intersection with the fractal set appeared 
Figure HI 

The key to the proof is a self-inducing structure with a scaling constant 
w^. We consider a region L' = uj~'^L and consider the first return map 

f(a;) =r™(^)(x) 

for X E L' where m{x) is the minimum positive integer such that T^'^^^x) G 
L' . For any x G L', the value m[x) = 1, 3 or 6. We can show that 

(1) uj'^f{uj-^x) = T{x) 

for X £ L. The proof is geometric, shown in Figure \E\ The return time 
m(x) = 3 in the open pentagon A' = a;~^A [this is marked A in the figure] 
and m{x) = 6 in the shaded pentagonal region D with three closed and two 
open edges. In the remaining isosceles triangle in uj^'^L (whose two equal 
edges are closed and the other open), the return time m(x) is 1. 

Note that the equation is valid for all x G L' . This makes the later 
discussion very simple. Unfortunately this is not the case for other quadratic 
values of 7 and we have to study the behavior of the boundary independently, 
see [3]- 

Let U be the 1-st hitting map to L' for x G L, i.e., U{x) = T"^^^\x) for 
the minimum non-negative integer m(x) such that T'^^^\x) £ L' . Note that 
f7 is a partial function, i.e., U{x) is not defined when there is no positive 
integer m such that T'"(x) G L' . Since 



T(x) 



x/C, X G A 

(x-l)/C xGr(Z)\A' 
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Figures. Self Inducing structure 
it is easy to make the map U explicit: 



U{x) 



f 
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L' 


(x-l)/C 
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T^{D) 




X 


G 


T\D) 
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T^{D) 




X 


G 


T^{D) 




X 


G 


T{D) 


Not defined 

V 


X 


G 


PoUPi 



where Pq is the largest open pentagon and Pi and P2 = Pi/C ^-re two second 
largest closed pentagons in Figure [6l 




Figure 6. Period Pentagons 
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C 1 1 



Set 

n — J n 1 _ 
to use later. 

We introduce a crucial map 5 which is the composition of the 1-st hitting 
map \J and expansion by w^, i.e. S{x) = uj'^U{x). Denote by 7r(x) the period 
of T-orbits of x G L and put 7r(x) = oo if x is not periodic by T. (We easily 
see tt{x) = 5 in Pq and 7r(x) = 10 in Pi U P2 unless x is the centroid of the 
pentagon.) Then if Tr{x) and 7r{S{x)) are defined and finite, then we see 
that Tr{S{x)) < Tr{x) which is a consequence of Equation ([T|). Therefore if 
71" (x) is finite then we have a decreasing sequence 

7r(x) > niSix)) > iriS^ix)) > ... 

of positive integers. This shows that there exists a positive integer k such 
that S^{x) is not defined. In this case we say that S-orbit of x in finite. 
We easily see that if S'-orbit of x G L is finite, then clearly 7r(x) is finite 
by Equation ([1]). Thus we have a clear distinction: x G L is T-periodic if 
and only if its S'-orbit is finite. Assume that x G L n gives an infinite 
S'-orbit. When U{x) is defined, we have U{x) = (x — c?rn(x))/C"'^^^ with 
m(x) = {0, 1, 2, 3, 4, 5} and di £ Q for all x G L. Thus we have 



dm 



By the assumption S^{x) is defined for A; = 1, 2, . . . and stays in L. Consider 
the conjugate map 4> which sends C C^- As (j){uj) = —l/co, we have 



rrii 



^2k(^2 E,Li m, ^ ^2{k-i+l) ^2 E,'=. rn, 

with = i;^>((ij) G (piQ)- Put A = max{|(i^| : G Q}- Then we have 

|0(S'=(x))|<|<^(x)| + ^^ 

Thus we have S^(x), (/)(S'^(x)) and their complex conjugates are bounded 
by a constant which does not depend on k. This implies that the sequence 
(S^(x))fc must be eventually periodic. 

Summing up, for a point x in its S-orbit is finite or eventually 

periodic. When it is finite then its T-orbit is periodic and when its S-orbit 
is eventually periodic then T-orbit is aperiodic. Thus we have an algorithm 
for a; G Z[/3] n L to tell whether its T-orbit is periodic or not. Since 

(x))| <—^ + ——^, 
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for any positive e, the right hand side is bounded by 

A 



e + 



a;2-l 

for a sufficiently large k. This means that under the assumption that there 
is an infinite /S-orbit, the set 

xeZ[C]nL : |0(x)| <e + — ^ 

— 1 

contains x with 7r{x) = oo. Since this set is finite, it is equal to 

f A 

B = \ x£'L[C]r\L : \(j){x)\ < 



Cj2_ 1 



for a sufficiently small e. Since there are only finitely many candidates in 
B, we obtain an algorithm to check whether an element x S Z[(^] n L with 
tt(x) = oo exists. In fact, all elements in B gives a finite S'-expansion, we 
are done. 

The same algorithm applies to jjZlQ with a fixed positive integer M. 
In this way, we can also show that points in are periodic. We can 

find aperiodic orbits in For example, one can see that 1/3 has an 

aperiodic T-orbit because its 5-orbit: 

1 nP oj^C^ _ 2£^ a;~\-i 

3' 3 ' 3 ' 3 3 ' 3 ' 3 '"' 

satisfies S^{l/3) = S'^il/S). 

It is crucial in the above proof that the scaling constant of the self-inducing 
structure is a Pisot number. Scaling constants of piecewise isometries often 
become Pisot numbers, moreover algebraic units. We discuss these phenom- 
ena in Section [3l 

2. Coding of aperiodic T-orbits 

Denote by A the set of all T-aperiodic points in L. By the proof of the 
previous section, we have 

A = {x G L I S''{x) is defined for aU A; = 1, 2, . . . }. 

We also have ^(A) C A. This means that for xi € A, there is a mj G 
{0,1,2,3,4,5} and Xj G A such that oj'^C,~^-^{xi — dmj = a^i+i G A for 
i = 1, 2, . . . . We therefore have an expansion 

(2) Xi — dfYii H ^ I H 9~ ( H ^ I '^r?i4 ~l~ n~ • • ■ 

Conversely a sequence {?rai}j=i,2,... defines a single point of Y. Therefore A 
must be a subset of the attractor Y of the iterated function system (IFS): 



i=0 
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an approximation of which is depicted in Figure 7(a) At this point we can 




0.2 0.4 0.6 . £ 

(a) All digits 




0.2 0.4 0.6 . 

(b) do,d2,d3,d5 



Figure 7. Attractors containing A 



assert that 2-dimensional Lebesgue measure of aperiodic points in L must 
be zero, because ~ 6.854 • • • > 6. 

We notice that the digits in Q are not arbitrarily chosen because the 
image of S must be in T{Z). Thus the digits di and ^4 appears only at the 
beginning in the expression of Equation ([2]). Therefore it is more suitable 
to study A n T{Z). The attractor 

(3) r = (-Ly + do) u (gy + d,) u (^^y + d,) u 

is depicted in Figure |7(b)[ 

This iterated function system satisfies OSC by a pentagonal shape K with 
whose vertices are 



o,-C\C,-C^-'-C' 

as in Figure [8j We confirm that the pieces K, 



Am 

+ dm do not overlap. 
We consider the induced system of {L,M,i^,T) to T{Z). Denote by T 
the first return map on T{Z). Then the induced system (T{Z),T) is the 
domain exchange of two isosceles triangle A and B depicted in Figure [9l 
The triangle A has two closed edges of equal length and one open edge, 
while B has one closed edge and two open edges of the same length. The 
open regular pentagon Pq and the triangle B move together by T and can 
be merged into a single shape. 
We see 



(4) 



r(x) 




X £ A 

X G T{Z) \ A. 
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T(B) 



T(Pc 



T(A) 




Figure 9. Induced Rotation T on T{Z) 

Again we find self- inducing structure with the scahng constant w^: 
(5) w^f (w-^x) = f{x) 

for all X G T{Z). This can be seen in Figure [TOl with a = uj~'^A, /3 = uj^'^B 
and R = uj^^Pq. This induced dynamics (T{Z),T) is essential in describing 
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Figure 10. Self Inducing Structure of (r(Z),T) 



the set A. 

Readers may notice that we can find a self-inducing structure by smaller 
scaling constant cj in Figure [7(b) | by taking two connected pieces. However 
this choice of inducing region is not suitable because the self-inducing rela- 
tion (with flipping) is measure theoretically valid, but has different behavior 
on the boundary. 

Let us introduce two codings. First is the coding of T-orbits of a point x 
in L in two symbols {0, 1}: d(x) = (^(r"(2;))„ G {0, 1}'^ where 



ip{x) 



For e.g., the d(l/3) =20110101011010101101101101 . . . The second coding 
is defined by d(x) = (V^(f "(a;)))^ E {a,b}^ for x e T{Z) where 




V'(x) 



a X G A 

b x£ T{Z) \ A. 



For a point x in T{Z) we have two codings by {a,b} and by {0, 1}. From 
Equation (HD, the two codings are equivalent through the substitution a — 
01, 6 — >■ 1. For a given coding of T-orbit by {0, 1}, there is a unique way 
to retrieve the coding of T-orbit by {a, 6}, because the symbol must be 
followed by 1. For e.g., r(l/3) = -2C~V3 G T{Z) is coded in two ways as: 

d{—2(^^/3) = ab a a ab a a ab ab ab ab . . . 
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and 



d(-2C"73) = 01 1 01 01 01 1 01 01 01 1 01 1 01 1 01 1 



Hereafter we discuss the coding d. Observing the trajectory of the region 
(x>~^(A) and lo^'^(T{Z) \ A) by the first return map by the iteration of T to 
the region lo^'^T{Z), it is natural to introduce a substitution ctq: 



aaba, 



baba. 



on {a, b}* and we have 

d{u}^'^x) = ao{d{x)) 

for X G T{Z). More generahy, fohowing the analogy of the previous section, 
the first hitting map to the region uo^'^(T{Z)) provide us an expansion of 
a point X E T(Z) exactly in the same form as ^ with restricted digits 
{do, (^2,^3,^5}. One can confirm that 



(6) 



7X + d„ 



[ao{d{x)) ^ 
a e ao{d{T{x))) 
ba^aoidij^x))) 
aba®aQ{d{T^{x))) 



m 
m 
m 
m 




2 
3 
5 



where ® is the concatenation of letters. Defining conjugate substitutions by 



O"! = aa^a ^, 0"2 = baaia ^b ^ and (T3 
(To (a) = aaba, 
cTi(a) = aaab, 
0-2(0) = baaa, 
0-3(0) = 0600, 



abaa^a ^b ^a ^, 



I.e., 



Mb) 
a,{b) 

Mb) 
Mb) 



one may rewrite 



tx + d„ 



' M^{x)) 
ai(d(r(x))) 
aM{T^{x))) 

ta3(d(T3(x))) 



baba 
abab 
baba 
abab 

m 
m 
m 
m 




2 
3 
5. 



We say that an infinite word y in {o, 6}^ is an S-adic Umit of ai {i 
if there exist yi G {a, b}^ for i = 1,2, . . . such that 



0,1,2,3) 



Um o-mi o 0-^2 o 0-^,3 o • • • o ami ive)- 



with rrii £ {0,1,2,3}. Since each element x £ T{Z) D A has an infinite 
expansion ([2]) with digits {doj '^2, c?3) c^s}) we find Xi € uj^'^T{Z) such that 



dfx) 



lim CTmi o 0-^2 o 0-^3 o 



o a, 



mi 



(d(x,)). 



This shows that d{x) is an S'-adic limit of {i = 0,1, 2, 3). 

Note that from the definition ^ of o-j, for a given S'-adic limit y there is 
an algorithm to retrieve uniquely the sequence {(Jmji- Checking first four 
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letters of y, we know the first letter of yi and to determine mi we need first 
6 letters. We can iterate this process easily. 

Summing up, we embedded the set A n T{Z) into the attractor Y' of an 
IFS ([3]) and succeeded in characterizing the coding of T-orbits of points in 
this attractor as a set of S'-adic limits on {ctq, cJi, (72, as}. However recalling 
that points in closed pentagons Pi and P2 are T-periodic and Y' is a non- 
empty compact set, we see from Figure [7(b)] that A is a proper subset of 
Y'. 

We wish to characterize the set of aperiodic points in Y' and its coding 
through d. Recalling the discussion in the previous section, if x G T{Z) 
has periodic T-orbits if and only if there exists a positive integer k such 
that S^{x) G Pq U Pi U ^2- The equivalent statement in the induced system 
(T{Z),T) is that x £ T{Z) is T-periodic if and only if there exists a positive 
integer k such that S^{x) € Pq U Pi. Note that we have: 



T{x) 



C, ^{x -p) +p X £ Pq 
C'^{x-q) + q xGPi 



where P = \ + i "V^^^^+^v^) ^^.ggp^ ^ _ i\J ^^^ ) is the center of Pq (resp. Pi) 

and consequently T^{x) = x holds for x G Pq U Pi. If x G T{Z) and x is 
T-periodic, then there exist Xi £ T{Z) such that xi £ PqU Pi and 

with m, G {0, 2, 3, 5}. Thus the set of T-periodic points in T{Z) consists of 
all the pentagons of the form 



(7) dm\ ~\~ ^2 ^'^"12 ^2 ^'^'"■3 ~^ ^2 ■ ■ ■ yd^l ~^ (j2£ ^ 

with J = 0, 1 and rrij G {0, 2, 3, 5}. From the self-inducing structure ([ll), it is 
easy to see that if two points x, x' are in the same pentagon of above shape 
and none of them is the center, then they have exactly the same periods. 
Moreover two T-orbits keeps constant distance, i.e., T'^(x) — T"'(x') = C^{x — 
x') for some integer s. The period is completely determined in |3]. We have 
all the periodic orbits in T{Z) and therefore have a geometric description of 
aperiodic points: 

A n T{Z) = T{Z) \ {All pentagons of the form 

Subtraction of these pentagons from T{Z) is described by an algorithm. 
The initial set is Dq = T{Z) \ Pq with two open and three closed edges as 
in the left Figure [TTJ The interior Inn(T'o) gives another feasible open set 
to assure the open set condition of the IFS of ([3]) . Inductively we define the 
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decreasing sequence of sets 

(/-m \ 

for i = 0, 2, . . . . Then Di consists of 4* pieces congruent to uj~'^''Dq without 
overlapping. Note that since 

Di = Do\ (Pi U u-^Po U ^ 

Di is obtained by subtracting from Dq one closed and two open regular 
pentagons as in Figure [TTl To generate -Dj+i, each 4* pieces in Di are 




Figure 11. Pentagon Removal Algorithm 

subdivided into 4 sub-pieces by subtracting three small regular pentagons. 
Clearly all regular pentagons of the shape ([7]) are subtracted by this iteration 
and we obtain 

oo 

AnT{Z) = f]Di. 

i=0 

This observation allows us to symbolically characterize aperiodic points in 
Y'. First, every point x of Y' has an address dmidm2 • • • S {do,d2,d-i,d5}^ 
by the expansion ([2|). The address is unique but for countable exceptions. 
The exceptional points forms the set of cut points of Y' having the eventually 
periodic expansion: 

dod2idor ^ dsdsid.r 

dsidor - d^id^r 

d2d2{dor ~ d^dsid^r 
in the suffix of its address, which is understood by Figure [12] where Kmn = 

^{^K + dn)+dm. 

Note that if a point x in T{Z) is periodic, then there exists a non-negative 
integer k such that S^{x) G Pq U Pi. Moreover, if x ^ Y' n T{Z), then there 
exists a non-negative integer k such that S^{x) G <9(Pi), because it can not 
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K55 K52 K23 K20 




Figure 13. T-periodic expansions 



be an inner point of Pq or Pi. In other words, such x must be located in 
the open edge of one of pieces of Dj.. From Figure [TU one can construct 
the following Figure [13] which recognize points of two open edges in d{DQ). 
For construction, we introduce a new symbol set {R, L} (right and left) to 
distinguish which open edge of is into focus. 

To read the graph and obtain the previous sequences, ignore {R, L} and 
substitute {0,2,3,5} with {do, ^2, c?3, ^5}. A point x £ Y' is periodic (or 
in the open edge of Dq) if and only if a suffix of the address dmidm2 • • • £ 
{do, d2,ds, d^}^ is in Figure [T3l Note that the points with double addresses 
are on the open edge of some Di and consequently their suffixes are read in 
Figure [T3l Figure [12] helps this construction. For e.g., the right open edge 
of consists of the left open edge of i^sa and the right open edge of K^q, 
therefore we draw outgoing edges from 5R to 3L and OR. 
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0,2,3,5 




Figure 14. Biichi automaton for periodic points in Y' 

As a result, the set of addresses of the points in A n T{Z) are recognized 
by a Biichi automaton which is the complement of the Biichi automaton of 
FigureO Here the double bordered states in FigurefT^are final states. Each 
infinite word produced by the edge labels {do) ^^2, c^3, t^s} on this directed 
graph is accepted, because it visits infinitely many times the final states. 
We do not give here the exact shape of its complement. It is known that 
complementation of a Biichi automaton is much harder than the one of a 
finite automaton, because the subset construction does not work (c.f. p9l 



Now consider the topology of {a, h}^ induced from the metric defined by 
r ^ ^ xiX2...,y = 2/12/2 •• • G {a,&}^- Take a fixed point 
w = (u)j)j=o,i,2,... £ {flj^}^ with (Jo{w) = w. This is computed for e.g., 
by lim„(7Q(a). The shift map y is a continuous map from {a, 6}^ to it- 
self defined by V{{'Wi)) = (wi+i). Letting Xa-Q be the closure of the set 
{V^lw) I n = 0,1,...}, we can define the substitutive dynamical system 
(Xg-Q,!^) associated with ctq. Since ctq is primitive the set X^-o does not de- 
pend on the choice of the fixed point and (Xg-Q , V) is minimal and uniquely 
ergodic (see [9}). Let r be the invariant measure of {Xuq,V). On the other 
hand, for the attractor Y' there is the self-similar measure v, i.e., a unique 
probability measure (c.f. Hutchinson [13]) satisfying 



for z^- measurable sets By/ in Y' . 

Theorem 3. The restriction ofT to Y' is measure preserving and (Y' , By / , z/, 
is isomorphic to the 2-adic odometer {Z,2,x i— )■ x + 1) as measure dynamical 



m)- 




me{0,2,3,5} 
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systems: 



I2 ^2 



(8) 



Y' 



-> Y' 



where </> : Z2 — t- y is almost one to one and measure preserving, which will 
be made explicit in the proof. Moreover the map 



p : X f-^ 



X — {x mod 4) 



from 7j2 to itself gives a commutative diagram: 

p 



(9) 



Z2 



Y' 



-> Y'. 



The above theorem may be read that , By, z^, T) gives a one-sided 
variant of numeration system in the sense of Kamae [15] . 



Proof. First we confirm that T is measure preserving. Denote by [dm^ , dm2 ; • • • ; 
the cyHnder set: 



(10) dra, + 



c 



^m2 



",2 



2 y dm2 ~l~ ^2 ( '^'"■3 ~^ n • • • I dfYii ~l~ o/? ^ I • • • 



c 



mil 



"mi, ^ J2l 



m\ ) d'ra2 ? 



. . . , dmi ] ) — 4 

[d3],f-l([(io]) = 

0,2,3, then iy{T-^{[dm^,dm2, ■ ■ ■,dme])) 



From Figure [TU[ we see 
[^2] but f-^{[d5]) inter- 



By the OSC, we have v^d. 
that f-i([d3]) = [ds], r-n[f^2]) 
sects both A and -B. Hence if mi 
4~^. By using the self-inducing structure in FigurefTOl we also have T~^{[d5d3]) 
[dod^], f-\[d5d2]) = [dods] and f-^id^do]) = [^0^2]. Thus if m2 = 0,2,3, 



then iy{T ^{[d^, d. 



m2 1 ■ ■ ■ 1 "'mi 



dmi])) = 4 . Repeating this, we can show that 



i^{T {[dmi,d 



dmi])) 



'7712 ) ■ ■ ■ ; "mi 

holds for all rui G {0, 2, 3, 5} but a single exception mi = m2 = ■ ■ ■ = mi = 5. 
Since i is arbitrary chosen, a simple approximation argument shows that T 
is measure preserving and (1"', By/, v, T) forms a measure dynamical system. 

Let us define a map r] from X^-q to Y' . Take an element z = X1X2 • • • S 
Xo-Q. Then each prefix X1X2 . . .Xi with £ > 3 is a subword of the fix point 
V of (To starting with a. Therefore there is a word y G {X,a,ba,aba} and 
zi G Xf^Q such that x = yiO'o{zi). It is easy to see from ([6]) that this y and 
zi are unique. 
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Iterating this we have Zi = yi^ia{zi+i) with yi G {X,a,ba,aba}, Zj+i S 
Xo-p and zq = z. Thus we have for any £, 

z = yio-o(y2 do (2/30-0 . . . yeic^oize)))) 

= yio-o(y2)o-o(y3) • • • (^o'^{ye)o-o{ze)- 

Define a map from {A, a, ba, aba} to Z by 

k{X) = 0, K{a) = 1, K{ba) = 2, K,{aba) = 3. 

Then Zi = yi+ia{zi+i) is equivalent to = o"K{yj+i)(^;j+i) and z is repre- 
sented as an S-adic limit: 

z = lim a^^iy^) o a^,^y^^ o • • • o cr^/j^ n(z^). 

for £ = 1,2,.... This gives a multiplicative coding d' : X^j^ — )• {ao, o"i, o"2) o"3} 
Let A' be the points of Xo-q whose multiplicative coding does not end up 
in an infinite word produced by reading the vertex labels of Figure [T5j Let 
us associate to z a 2-adic integer i{z) = — X^^^q ^iVi)'^'^^ ^ ^2- The map t 
is clearly bijective bi-continuous and the value i{z) is also called the multi- 
plicative coding of z. We write down first several iterates of V on the fix 
point of do, to illustrate the situation: 



(ToCToCroCro • • • 


A -0000 


a^a^a^a^ . . . 


4 -3333 


0'2CT3Cr3Cr3 . . . 


A -2333 


aia^a^as . . . 


4 -1333 


aoa^asas . . . 


A -0333 


cr3CJ2Cr3Cr3 . . . 


A -3233 



One can see that the following commutative diagram (jlip holds. 

Y ^ , X 

^ CTQ ' ^ (TO 

(11) 
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Therefore (Xq-q, V) is topologically conjugate to the 2-adic odometer (^2,2; 1— )• 
X + 1). Here the consecutive digits {0, 1} in Z2 are glued together to give 
{0, 1, 2, 3} = {0, 1} + 2{0, 1}. Indeed, do satisfies the coincidence condition 
of height one in the sense of Dekking |25|f9j and above conjugacy is a conse- 
quence of this. (Z2, X I— )• x + 1) is a translation of a compact group Z2 which 
is minimal and uniquely ergodic with the Haar measure of Z2. Moreover one 
can confirm that l preserves the measure and (X^q, V) and (Z2, x 1— )• x + 1) 
are isomorphic through l as measure dynamical systems. In view of ([6]), we 
define 



'0 


i 


= X 


2 

< 


i 


= a 


3 


i 


= ha 


.5 


i 


= aha 



and the map rj 
(12) = 




^€(2/2) + 



+ 



Then rj is clearly surjective, continuous, and measurable because both r 
and V are Borel probability measures. Since the set of points with double 
addresses is on the open edge, the map 77 is bijective from A' to A n T[Z). 
Since d(T(x)) = V{d{x)), we have a commutative diagram: 

V 



(13) 



A' 



Anr(^) 



A' 



AnT(z). 



From Figure [T^ it is easy to see that the set V of T-periodic points in Y' 
is measure zero by i.e., v{A n T{Z)) = v{Y' n T{Z)) = 1, because the 
number of words of length n in Figure [T3] is 0(2"). Similarly as the Perron- 
Frobenius root of the substitution (Tq is 4 and the number of words of lengths 
n in Figure US] are 0(2"), we see that r(A') = r(X^o) = 1. From ([13]) the 



pull back measure i' o r] ^ of Xa-^ is invariant by V, we have r 



u o 



7] 1 by 



unique ergodicity. Therefore by taking (/> = ry o t, we have the commutative 
diagram ([8]) with measure zero exceptions. Let V be a map from Xq-q to 
itself which acts as the shift operator on the multiplicative coding d', i.e.. 



(14) 



n2'~'n-i ■ ■ ■ 



for d'{z) 
A' 



(^ni • • • 



Then we see that 



V 



A' 



AnT(Z) 



^ AnT(z). 



and the commutative diagram ([9]) is valid but for measure zero exceptions. 

□ 
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Corollary 4. Each aperiodic point x G A n T{Z), the T-orbit of x is uni- 
formly distributed in Y' with respect to the self similar measure v. 

Proof. In the proof of Theorem [3] the map r/ is bijective form A' to AnT(2). 
Therefore if x € A n T{Z), then there exists a unique element in z £ X^g 
with r]{z) = X. Therefore there exist an element zq G Z2 such that (j){zo) = x. 
The Haar measure /Li2 on Z2 is given by the values on the semi-algebra: 

Ai2([co,ci, . . . ,Q_i]) = 4"^ 

for each cylinder set [cq, ci, . . . , = {2/eZ2 I y = Etdci4^ (mod 4^)}. 
Since (Z2, x 1— t- x + 1) is uniquely ergodic, the assertion follows immediately 
from the commutative diagram ([8]). □ 
Not all points in Y' gives a dense orbit as we already mentioned that 
Any' is a proper dense subset of Y' . There are many periodic points in Y' 
as well. This gives a good contrast to usual minimal topological dynamics 
given by a continuous map acting on a compact metrizable space. 

Corollary 5. Each aperiodic point x G A, the T-orbit of x is dense in the 
set X. 

Proof. It is clear from the fact that (T(Z),T) is the induced system of 

(L,r). □ 

One can construct a dual expansion of the non-invertible dynamics (Y' , S) 
by the conjugate map (/> : C — ^ iii Gal(Q(C)/Q) and then make a natural 
extension: an invertible dynamics which contains {Y',S). The idea comes 
from symbolic dynamics. We wish to construct the reverse expansion of (I12p 
to the other direction. To this matter, we compute in the following way: 



^C(s/2) 1 = d^iy-i) + ( d^{y4) + 




Therefore it is natural to introduce a left 'expansion': 



up' ( ( tjp 



with ifc G {0,2,3,5}. As this expression does not converge, we take the 
image of because 4>{oS) = —l/u. Let us denote by tij^, = <t){di^). Then the 
expansion 

A-2ii //•-2i2 /A-2j3 



\ I ^;2 



converges and the closure of the set of such expansions gives a compact set 
y depicted in figure [I6l 
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Figure 16. The dual attractor y 



Of course the set is an attractor of the IFS: 



y 



1 



c 



{y -uo)u^{y- U2) u^{y-u3)u -^{y - «5). 



Combining y we can construct a natural extension of (Y' , S) as: 



Y' xy3 {ri, 9) A 



{ii-di)u' c-'%e - <p{di)) 







eY' xy 



On the other hand (7j2^ p) have a natural extension: 



, N , ^ P f X — ix mod 4) y + (x mod 4) , „ , x 

Z2X [0,l)9(x,y)4 ^ -] eZ2x[0,l) 



and two systems are isomorphic both as topological and measure theoretical 
dynamics: 



(15) 



where 0' is given as: 



Z2 X [0, 1) 

Y' xy 



Z2 X [0, 1) 

0x</>' 

Y' xy. 



^Xi4 ' ^ gxi{gx2{9x3{- ■ ■))) 



1=1 



where gQ{x) = {x - uq) / uj"^ , gi{x) = {x - U2)C, / uP' , g2{x) = {x - usX 
and gsix) = {x - u^)/^:'^. 

From this 'algebraic' natural extension construction, we can characterize 
purely S'-periodic points in Y' fl Q(C)- 
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Theorem 6. A point y in Y'nQ{C) has purely periodic multiplicative coding 
with four digits do, o"2, crs, us if and only if {y, 0(y)) £ Y' x y. 

This is an analogy of the results [T4j for /3-expansion. The proof below is 
on the same line. 

Proof. As uj is an algebraic unit and di G the denominator of gi{y) is 
the same as that of y for i = 0, 1, 2, 3. Therefore the module y £ Ai = jjZ[C,] 
is stable by gi for some positive integer M. Note that points y £ M with 
{y, (j){y)) G y X 3^ is finite, because y, (j){y) and their complex conjugates are 
bounded in C. One can confirm that the map S becomes surjective from Ai 
to itself. For a finite set, surjectivity implies bijectivity. Therefore a point 
y £ M. with {y,4>{y)) £ Y' x y produces a purely periodic orbit. On the 
other hand if x has purely periodic multiplicative coding, it is easy to see 
iy,cP{y))GY' xy. □ 



3. Other self-similar systems 

Pisot scaling constants appear in several important dynamics. For ir- 
rational rotations (2IET), it is well known that scaling constants of self- 
inducing systems must be quadratic Pisot units. A typical example Figure [2] 
was shown in the introduction. They are computed by the continued fraction 
algorithm as fundamental units of quadratic number fields. Poggiaspalla- 
Lowenstein-Vivald [24j showed that the scaling constant must be an alge- 
braic unit for self-inducing uniquely ergodic lET. When the scaling constant 
of self-inducing lET is a cubic Pisot unit, we have further nice properties 

[a da [20]. 

A necessary condition that 1-dimensional substitutive point sets give point 
diffraction is that the scaling constant is a Pisot number [7]. Suspension 
tiling dynamics of such substitution is conjectured to have pure discrete 
spectrum if the characteristic polynomial of its substitution matrix is irre- 
ducible. For higher dimensional tiling dynamics the Pisot (or Pisot fam- 
ily) property is essential to have relatively dense point spectra, see for e.g. 

Pisot scaling properties seem to extend to the case of piecewise isometrics. 
To conclude we present some examples, though we do not make a systematic 
study. 

It is already observed in [16^ [3] that Pisot scaling constants appear in our 
problem if 9 is the n-th root of unity for n = 4, 6, 8, 10, 12 in the same way 
as we did in ra = 5 but in a more involved manner. In each case they are 
quadratic Pisot units. What about if A = — 2 cos(0) is cubic? In this case, the 
dynamics of Conjecture [T] are embedded into the piecewise affine mapping 
acting on (R/Z)^ which is harder to visualize. Instead let us consider formal 
analogies of piecewise isometrics generated by cubic n-th fold rotation in the 
plane. At the expense of losing connection to Conjecture [U we find many 
Pisot unit scaling constants! Being an algebraic unit is natural and may 
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be explained from invertibility of dynamics. However we have no idea why 
the Pisot numbers turn up or even how to formulate these phenomena as a 
suitable conjecture. 

3.1. Seven-fold. We start with 7-fold case. Both pieces are rotated clock- 
wise by 47r/7 as in Figure [T71 The triangle is rotated around A and the 
trapezium around B. The first return map to a region and a smaller region 
with the same first return map (up to scaling) are described. Unlike the 
five fold case, returning to the subregion does not cover the full region. A 
simple consequence is that there are infinitely many possible orbit closures 
for non-periodic orbits in the system. The scaling constant a w 5.04892 is 
a Pisot number whose minimal polynomial is x — + 5x - 1. Fi gure [m 
shows how this remaining space can be filled in. As this region is already a 
little small we will zoom in and now consider just this induced sub-system in 
Figure [T9l The smaller substitutions are easier to see as there are two scal- 
ings giving the same dynamics (A and B). The scaling constant f3 w 16.3937 
for these subregions is the Pisot number associated to — 17x^ -|- lOx — 1. 
The proof that the remaining substitutions work is shown in Figure [20l The 
first return map to the two lower triangles is shown. The same dynamics 
occur on a smaller region. The orbit of the smaller region covers all the 
regions left out of Figure [19] and so the substitution rule from that figure 
is now complete. The scaling constant for this triangle is a. This gives an 
example of recursive tiling structure by Lowenstein-Kouptsov- Vivaldi [18] . 
Knowing that every aperiodic orbits are in one of the above self-inducing 
structures, we can show that 

Theorem 7. Almost all points of this 7-fold lozenge have periodic orbits. 



The argument is similar to that given around Figure 7(a) We easily find 
decreasing series X„ of union of polygons satisfying /i(aX„_|_i) < a^/i(X„) 
(or fi{/3Xn+i) < /3^/i(X„)) which cover all self-inducing structures. 

The fundamental units of the maximal real subfield Q(cos(47r/7)) of the 
cyclotomic field Q(C7) are given by b and b — 1 where b = 1/(2 cos (37r/7)) ~ 
2.24698. Here b is the Pisot number satisfying x^ — 2x'^ — x + 1. We see 
that a = b'^ and /3 = b^/{b — 1)^ and thus a and /3 generates a subgroup 
of fundamental units of Q(cos(47r/7)). Note that both ^/a = b and = 
6^/ (b—l) are Pisot numbers but 6—1 is not. Our piecewise isometry somehow 
selects Pisot units out of the unit group! 

3.2. Nine-fold. The next example is 9-fold case in Figure [211 Both pieces 
are rotated anti-clockwise by 47r/9, the triangle around A and the trapezium 
around B. The first return map (A) to the triangle is also shown. In addition 
the same dynamics are found on a smaller piece of the map. Like the 7- fold 
shown in Figure [17] this does give a full description of the dynamics, but it 
is not A. The scaling constant 7 ~ 8.29086 is a Pisot unit defined by 
x'^ — 9x^ + 6x — 1. Unfortunately in this case we were not able to find a 
complete description of the scaling structure. 
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First return to region 




SeLf-inducing map 




Figure 17. A seven- fold piecewise isometry. 



The fundamental units of Q(cos(47r/9)) are b and 6^ — 26 — 1 where b = 
1/(2 cos(47r/9)) w 2.87939 is a Pisot number given by — + 1. We have 
7 = 6^ and are expecting to find another Pisot unit 6^/(6^ — 26—1) ~ 5.41147 
(or its square) as a scaling constant in this dynamics, which would give an 
analogy to the seven-fold case. 
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Figure 21. A nine fold piecewise isometry. 
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